AIC, BIC, Bayesian evidence and a notion on simplicity of cosmological model 



O 
(N 



Marek Szydlowski* 

Department of Theoretical Physics, Catholic University of Lublvn, Al. Raclawickie I4, 20-950 Lublin, Poland and 
Mark Kac Complex Systems Research Centre, Jagiellonian University, Reymonta 4, 30-059 Krakow, Poland 

Aleksandra Kurek^ 

Astronomical Observatory, Jagiellonian University, Orla 171, 30-244 Krakow, Poland 

Recent astronomical observations indicate that the Universe is in the phase of accelerated expan- 
sion. There are many cosmological models which explain this phenomenon, but should we prefer 
those models over the simplest one - ACDM model? According to the Occam's razor principle if 
all models describe the observations equally well we should prefer the simplest one. We consider 
the model comparison methods which involve such rules: the Akaike information criterion (AIC), 

00 ' Bayesian information criterion (BIC) and Bayesian evidence to compare the ACDM model with its 
' generalisation where the interaction between dark matter and dark energy is allowed. The analyses 

based on the AIC and Bayesian evidence indicate that there is only a weak evidence in favour of 
the ACDM model over its generalisation, while those based on BIC quantity indicate the strong 
' evidence in favour the simpler model. We also calculate some quantity which measure the effective 

, number of model parameters that the given data can constrain. This value is used to compare the 

concordance LCDM model with its generalization basing on the extended interpretation of continu- 
ity condition-interacting ACDM cosmology. We conclude that data set are not enough informative 
to constrain all parameters allows to vary in the models. 

2 . I. INTRODUCTION 

1 ^ I ' Recent observations of type la supernova (SNIa) provide the main evidence that current Universe is in an accelerating 
phase of expansion [1]. Cosmic microwave background (CMB) data indicate that the present Universe has also 
negligible space curvature [2]. Therefore if we assume the Friedmann-Robertson- Walker (FRW) model in which 

^ [ effects of nonhomogeneities are neglected, then the acceleration must be driven by dark energy component X (matter 
fluid violating the strong energy condition px + 3px > 0). This kind of energy represents roughly 70% of the matter 
content of the current Universe. Because the nature as well as mechanism of the cosmological origin of the dark energy 
component are unknown some alternative theories try to eliminate the dark energy by modifying the theory of gravity 
itself. The main prototype of this kind of model is covariant brane models based on the Dvali-Gabadadze-Porrati 
(DGP) model [3] as generalized to cosmology by Deffayet [4]. The simplest explanation of a dark energy component is 
the cosmological constant with effective equation of state p — —p but appears the problem of its smallness and hence 
I its relatively recent dominance. Although the ACDM model offers possibility of explanation of observational data it 
• • . is only effective theory which contain the enigmatic theoretical term - the cosmological constant A. Other numerous 
candidates for dark energy description have also been proposed like to evolving scalar field [5] usually referred as 
quintessence, the phantom energy [6, 7], the Chaplygin gas [8] etc. Some authors believed that the dark energy 
' problem belongs to the quantum gravity domain [9] . 
5^ ', These theoretical models are consistent with the observations, they are able to explain the phenomenon of the 
accelerated expansion of the Universe. But should we really prefer such models over the ACDM one? To answer this 
question we should use some model comparison methods to confront existing cosmological models having observations 
at hand. 

Let us assume that we have N pairs of measurements (j/i, Xi) and that we want to find the relation between the y 
and X quantities. Suppose that we can postulate k possible relations y = fi{x,0), where 9 is the vector of unknown 
model parameters and i — 1, . . . ,/c. With the assumption that our observations come with uncorrelated gaussian 
errors with mean fii = and standard deviation at the goodness of fit for the theoretical model is measured by the 
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quantity given by 

^ ^ = -21nL, (1) 

^=l 

where L is the hkeUhood function. For the particular family of models /; the best one minimize the quantity, which 
we denote 0). The best model from our set of fc models ]\(x, 9), . . . , fk{x, 0) could be the one with the smallest 
value of quantity. But this method could give us misleading results. Generally speaking for more complex model 
the value of is smaller, thus the most complex one will be choose as the best from our set under consideration. 

A clue is given by so called Occam's razor principle: "If two models describe the observations equally well, choose 
the simplest one." This principle has aesthetic as well as empirical justification. Let us quote simple example which 
illustrate this rule [10] and present it in Figure 1. We see the black box and the white one behind it. One can 
postulate two models: 1. There is one box behind the black box, 2. There are two boxes of identical height and colour 
behind the black box. Both models explain our observations equally well. According to the Occam's razor principle 
we should accept the explanation which is simpler so that there is only one white box behind the black one. Is not it 
more probable that there is only one box than two boxes with the same height and colour? 



OBSERVATION 




HYPOTHESIS 1 HYPOTHESIS 2 




FIG. 1: Illustration to the example explaining the Occam's razor principle. 

We could not use this principle directly because the situations when two models explain the observations equally 
well are rare. But in the information theory as well as in the bayesian theory there are methods for model comparison 
which include such rule. 

In the information theory there are no true models. There is only reality which can be approximated by models, 
which depend on some number of parameters. The best one from the set under consideration should be the best 
approximation to the truth. The information lost when truth is approximated by model under consideration is 
measured by so called KuUback-Leibler (KL) information so the best one should minimizes this quantity. It is 
impossible to compute the KL information directly because it depends on truth which is unknown. Akaike [11] found 
approximation to the KL quantity which is called the Akaike information criterion (AIC) and is given by 

AIC = -21n£ + 2d, (2) 

where C is the maximum of the likelihood function and d is the number of model parameters. Model which is the 
best approximation to the truth from the set under consideration has the smallest value of the AIC quantity. It is 
convenient to evaluate the differences between the AIC quantities computed for the rest of models from our set and 
the AIC for the best one. Those differences (Aaic) are easy to interpret and allow a quick 'strength of evidence' for 
considered model with respect to the best one. The models with < Aaic < 2 have substantial support (evidence), 
those where 4 < Aaic ^ 7 have considerably less support, while models having Aaic > 10 have essentially no support 
with respect to the best model. 

It is worth noting that the complexity of the model is interpreted here as the number of its free parameters that can 
be adjusted to fit the model to the observations. If models under consideration fit the data equally well according to 
the Akaike rule the best one is with the smallest number of model parameters (the simplest one in such an approach) . 
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In the Baycsian framework the best model (from the model set under consideration) is that which has the largest 
value of probability in the light of data (so called posterior probability) [12] 

where P(Mj) is a prior probability for the model Mj, P{D) is normalization constant [P{D) = Yli=i -P(-0|Mj)P(Mi)], 
D denotes data. P{D\Mi) is the marginal likelihood, also called evidence 

P{D\Mi) = J P{D\9,Mi)P{e\Mi)de = Ei, (4) 

where P(D\$,Mi) is likelihood under model i, P{6\Mi) is prior probability for under model i. 

Let us note that we can include the Occam's razor principle by assuming the greater prior probability for simpler 
model, but this is not necessary and rarely used in practice. Usually one assume that there is no evidence to favor 
one model over another which cause to equal value of prior for all models under consideration. It is convenient to 
evaluate the posterior ratio for models under consideration which in the case with flat prior for models is reduced to 
the evidence ratio (so called the Bayes factor B). The interpretation of the logarithm of Bayes factor values are as 
follows: < InS < 2 as a weak, 2 < InB < 6 as a positive, 6 < InB < 10 as a strong and In i? > 10 as a very 
strong evidence in favor better model. This quantity involves the Occam's razor rule. Let us simplify the problem to 
illustrate how this principle works here [10, 13]. 

Assume that P(6\D, M) is the non normalized posterior probability for the vector 9 of model parameters. In this 
notation E = J P{0\D,M)d0. Suppose that posterior has a strong pick in the maximum: ^mod- It is reasonable to 
approximate the logarithm of the posterior by its taylor expansion in the neighbour of ^mod so we finished with the 
expression 

P{e\D, M) = P{eMOB\D, M) exp [-{6 - eMOufC-\0 - ^mod)] , (5) 



where [C ^] . . = — ^ '"o^^'ae"'^^"' - - • Posterior is approximated by the gaussian distribution with 

the covariance matrix C and the mean ^mod- The expression for the evidence take a form E = 
P{0mod\D,M) J exp [— (^ — ^mod)^C~^(^ — ^mod)] d 0. Because the posterior has a strong pick near the maxi- 
mum, the most contribution to the integral comes from the neighbour close to ^mod- Contribution from the other 
other region of can be ignored, so we can expand the limit of the integral to whole i?**. With this assumptions one 
can obtain E = (27r) ^ Vdet CP(^mod \D,M) = (27r)7 v^detCP(£)|^MOD, M)P(^mod|M). Suppose that the likelihood 
function has sharp pick in 9 and the prior for is nearly flat in the neighbour of 0. In this case 9 = ^mod and the 

expression for the evidence takes the form E = £(27r)^-\/det CP(^|M). The (27r)5\/det CP(^|M) quantity is called 
the Occam factor (OF). When we consider the case with one model parameter with flat prior P{9\M) = the 
Occam factor 0F= which can be interpreted as the ratio of the volume occupied by the posterior to the volume 
occupied by prior in the parameter space. The more parameter space wasted by the prior the smaller value of the 
evidence. It is worth noting that the evidence does not penalize parameters which are unconstrained by the data [14]. 

As the evidence is hard to evaluation an approximation to this quantity was proposed by Schwarz [15] so called 
Bayesian information criterion (BIG) and is given by 

BIC = -21n£ + 2(ilnA/-, (6) 

where iV is the number of the data points. The best model from the set under consideration is this which minimize 
the BIG quantity. It is also convenient to analyse the difference between BIG quantities for the rest models from our 
set with the BIG for the best one. These differences could be interpreted in the following way: < Abic < 2 as a 
weak, 2 < Abic < 6 as a positive, 6 < Abic < 10 as a strong and Abic > 10 as a very strong evidence in favour of a 
better model. 

One can notice the similarity between the AIC and BIG quantities though they come from different approaches to 
model selection problem. The dissimilarity is seen in the so called penalty term: ad, which penalize more complex 
models (complexity is identifled here as the number of free model parameters). One can evaluated the factor by which 
the additional parameter must improve the goodness of fit to be included in the model. This factor must be greater 
than a so equal to 2 in the AIG case and equal to In N in the BIG case. Notice that the latter depends on the number 
of the data points. 

It should be pointed out that presented model selection methods are widely used in context of cosmological model 
comparison [13, 14, 16-35]. 
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We should keep in mind that conclusions based on such quantities depend on the data at hand. Let us mention 
again the example with the black box. Suppose that we made a few steps toward this box that we can see the 
difference between the height of the loft and right side of the white box. Our conclusion change now. 

Let us quote example taking from [19]. Assume that we want to compare the Newtonian and Einsteinian theories 
in the light of the data coming from laboratory experiment where general relativistic effects are negligible. In this 
situation Bayes factor between Newtonian and Einsteinian theories will be close to unity. Whereas comparing the 
general relativistic and Newtonian explanations of the deflection of a light ray that just grazes the sun's surface give 
the Bayes factor ~ 10^° in the favor of the first one (and even greater with more accurate data). 

Having this in mind an interesting supplement to the above considerations seems to be quantity which measure the 
effective number of model parameters (C(,) that the given data set can constrain [26, 27]. It can be computed using 
the relation 

C,=^)-xHO), (7) 

where the mean is taken over the posterior pdf for 9 and 9 is the mode of the posterior pdf (different choices are 
also possible). As have been shown this quantity correspond to the number of parameters for which the width of the 
posterior probability distribution is significantly narrower than the width of the prior probability distribution [26]. 
These parameters can be considered to have been well measured by the data given our prior assumption in the model. 
It helps to determine if the data is informative enough to measure the parameters under consideration. 

We share with George Efstathiou opinion [36-38] that there is no sound theoretical basis for considering the dy- 
namical dark energy, where as we are beginning to see an explanation for a small cosmological constant emerging from 
more fundamental theory. In our opinion the ACDM model has the status of satisfactory effective theory. Estathiou 
argued why the cosmological constant should be given higher weight as a candidate for dark energy description than 
dynamical dark energy. In this argumentation Occam's razor is used to point out a more economical model explaining 
the observational data. On the other hand Bicsiada advocated the use of Akaike information criterion which favour 
rather a dynamical model of dark energy (quintessence model) [20]. In our opinion quantification of Occam's razor 
by computing Bayesian evidence and Bayesian information criterion is more suitable as the AIC is unadequate for 
many statistical problems [39] . 

The main aim of this paper is to compare the simplest cosmological model - the ACDM model - with its gen- 
eralisation where the interaction between dark energy and dark matter sector is allowed using methods described 
above. 



II. INTERACTING MODEL 



The interacting interpretation of continuity condition (conservation condition) was investigated in the context of 
the coincidence problem since the paper Zimdahl [40], for recent developments in this area see Olivares et al. [41, 42], 
see also Le Delliou et al. [43] for discussion recent observational constraints. 

Let us consider two basic equations which determine the evolution of FRW cosmological models 

l = -lip-^^P) (8) 



p=-3H{p + p). (9) 

Equation (8) is called the accelerated equation and equation (9) is the conservation (or adiabatic) condition. Equation 
(8) can be rewritten to the form analogous to the Newtonian equation of motion 

dV 

where V = V{a) is potential function of the scale factor a. To evaluate V{a) from (10) via integration by parts it is 
useful to rewrite (9) to the new equivalent form 

|(^=')+4(«=) = o. (U) 

Prom (8) we obtain 

|^ = ^(, + 3p)d(a^). (12) 
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It is convenient to calculate pressure p from (11) and then substitute to (12). After simple calculations we obtain 
from (12) 



dV 
da 



a'^+pd{a') 
da 



(13) 



Therefore 



V = - 



pa 



(14) 



In formula (14) p means effective energy density of the fluid filled the Universe. 

Wc find very simple interpretation of (8): the evolution of the Universe is cqiiivalent to motion of the particle 
of unit mass in the potential well parameterized by the scale factor. In the procedure of reduction of the problem 
of FRW evolution to the problem of investigation dynamical system of a Newtonian type we only assume that the 
effective energy density satisfies the conservation condition. We do not assume the conservation condition for each 
energy component (or noninteracting matter sectors). 

Equations (8) and (9) admit the first integral which is usually called the Fricdmann first integral. This first integral 
has a simple interpretation in the particle-like description of the FRW cosmology, namely energy conservation 



- + V{a)=E = - 

where k is the curvature constant and V is given by formula (14). 
Let us consider the universe filled with two components fluid 



k 



P = Pni+PX, P = + WxpX, 



(15) 



(16) 



where Pm means energy density of usual dust matter and px denotes energy density of dark energy satisfying the 
equation of state px = wxpx, where wx = Wx{a). Then equation (11) can be separated on dark matter and dark 
energy sectors which in general can interacts 



|(,.a3) + 0*|(a3) 
j^ipxa^)+wxia)pxj^ia^)- 



r 
-r 



In our previous paper [44] it was assumed that 



which able us to integrate (17) which gives 



a 



C a 
a'^ n 



(17) 



(18) 



^ + ^{l+wx{a))px=-aa^-^. 

da a 

The solution of homogeneous equation (19) can be written in terms of average Wx{a) as 

-3(l+wT{a)) 



px = PXfiO. 



where 



Wx (a) 



j Wx{a)d{ln a) 
d{ln a) 



The solution of nonhomogeneous equation (19) is 



px 



-3(l+«7x(a)) 



Cx 



a 



3(i+w3r(o)) ■ 



(19) 



(20) 



(21) 



(22) 
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Finally we obtain 



Peff = + S^=Pm+Px= (23) 



a3 n + 1 a3(i+iBT(a)) 



^ an-l+3«.x(a)^^ 



3(l+t«3r(a)) , ^--'^ 



£j3(l+t;^(a)) ■ 



The second and last terms origin from interaction between dark matter and dark energy sectors. 

Let us consider the simplest case of wx{a) =const= wx{a)- Then integration of (22) can be performed and we 
obtain 

Cm Cx Gnt 

Z'c" = ^ + + ^ (24) 

where Cint = ^ ~ n-3wx ' ^^^^ '^^^^ ^® obtain one additional term in peff or in the Friedmann first integral 
scaling like a^~". It is convenient to rewrite the Friedmann first integral to the new form using dimensionless density 
parameter. Then we obtain 

= "-"'"^^ + + + + + + + zf^'+^-i (25) 

Note that this additional power law term related to interaction can be also interpreted as the Cardassian or 
polytropic term [45, 46] (one can easily show that the assumed form of interaction always generates a correction of 
type a"^,m = 1 — n, in the potential of the ACDM model and vice versa). Another interpretation of this term can 
origin from Lambda decaying cosmology when Lambda term is parametrized by the scale factor [47]. 

In the next section we draw a comparison between the above model with the assumption that TUx{a) = const = — 1 
and the ACDM model. 



III. RESULTS 



To compare the ACDM model with the interacting ACDM model we use the SNIa data, constraints from the CMB 
shift parameter, constraints from the SDSS parameter A as well as H{z) observational data. 
We use A^i = 192 SNIa data [48-50]. In this case the likelihood function has the following form 



LsN oc exp 



' Ni 



the 



obs\2 
i ) 



(26) 



where CTj is known, jif'^ = mi~ M (mj-apparent magnitude, M-absolute magnitude of SNIa), f//^^"^ = 5 log^Q DLi + 
M, M. = — 51ogiQ Hq + 25 and Dn = HodLi, where rft, is the luminosity distance, which with the assumption k = 
is given by cIm = (1 + Zt)cJ^" jj^. 

We also include information obtained from the CMB data. Here the likelihood function has the following form 



Lr oc exp 

where R is the so called shift parameter, R^^<^°'^ 



{K 



the 



^obs)2 



L70±0.03 for z^e 



(27) 

1089 [2, 51]. 



/o'" ^dz, and 

As the third observational data we use the measurement of the baryon acoustic oscillations (BAO) from the SDSS 
luminous red galaxies [52]. In this case the likelihood function has the following form 



Lji (X exp 



^^theor _ ^obs^2 



(28) 



where A^^^°' = ^/iKfi 



(^) 
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and A°^^ = 0.469 ± 0.017 for za = 0.35. 



ZA Jo H{z)""^ 

Finally we used the observational H{z) data {N2 = 9) from [53] (see also [54, 55] and references therein). This data 
based on the differential ages ( J|) of the passively evolving galaxies which allow to estimate the relation H{z) = ^ = 
— ^ . Here the likelihood function has the following form 



Lh oc exp 



{H{z,) - H,{z,)f 

.i=l 
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where H{z) is the Hubble function. Hi, Zi are observational data. 
The final likelihood function is given by 

L = LsnLrLaLh, (29) 

with the number of data points N = 192 + 1 + 1+9. 

To obtain the values of AIC and BIC quantities we perform the = ~2 In L minimization procedure after marginal- 

ization over the Hq parameter in the range < 60,80 >. The values of Bayesian evidence as well as effective number 
of model parameters were obtained using the MCMC algorithm called nested sampling [56] which implementation to 
the cosmological case is available as a part of the CosmoMC code [57, 58] called CosmoNest [29, 30, 34, 59] which was 
changed for our purpose. We assume flat prior probabilities for the model parameters in the range: Hq G< 60, 80 >, 
f^m,o G< 0, 1 >, r^int G< —1,1 >, 3 — n = m G< —10, 10 >. The values of evidence and Cb were averaged from the 
eight chains. Results are presented in Table I. 

TABLE I: Values of AAIC, ABIC, In B with respect to the ACDM model and for both models. 



model 


AAIC ABIC InS Ct 


ACDM 
Int ACDM 


0.7 ± 0.1 
1.28 7.9 1.7 ± 0.4 1.5 ± 0.2 



Analysing the value of difference between the AIC quantity one can conclude that there is substantial support for 
the interacting ACDM model with respect to the ACDM model while the difference between the BIC quantity to the 
conclusion that there is strong evidence against the more complex model. The logarithm of the Bayes factor quantity 
confirm the first conclusion there is only weak evidence in favor of the ACDM model. Analysing the effective number 
of model parameters we can conclude that the data sets used in evaluation are not enough informative to constrain 
all model parameters: in the ACDM case there is nearly one parameter which is well constrained by the data while in 
the case of its generalisation the number of effective parameters is equal to 1.5 (whereas the number of all considered 
parameters is equal to 4 here). 



IV. CONCLUSION 



We presented the methods of model comparison coming from information as well as bayesian theory. Both of them 
include the Occam's razor principle which states that if two models describe the observations equally well we should 
choose the simpler one. According to Akaike and Scliwarz rule the model complexity is interpreted in the term of 
number of free model parameter while according to the Bayesian evidence more complex model wast greater volume 
of the parameter space. Finally we present the quantity which measures the effective number of model parameters 
which the data used in analysis can constrain. This quantity is also called the Bayesian complexity and reduce to 
the number of free model parameters when the data set is highly informative. We use those methods to answer 
the question if we should prefer the generalisation of the simplest ACDM model where the interaction between dark 
matter and dark energy sector is allowed. The AIC and bayesian evidence give similar results: there is only weak 
evidence to favor the ACDM model over the interacting ACDM one. This is with contrary with the conclusion from 
the analysis of the BIC quantity: here there is a strong evidence to favor the ACDM model. The analysis of the 
Bayesian complexity gives that data sets considered carry not enough information to constraint all parameters allows 
to vary in the models considered. 



Acknowledgments 

This work has been supported in part by the Marie Curie Actions Transfer of Knowledge project COCOS (contract 
MTKD-CT-2004-517186). 



[1] A. G. Riess et al. (Supernova Search Team), Astron. J. 116, 1009 (1998), astro-ph/9805201. 
[2] D. N. Spergel et al. (WMAP) (2006), astro-ph/0603449. 

[3] G. R. Dvali, G. Gabadadze, and M. Porrati, Phys. Lett. B485, 208 (2000), hep-th/0005016. 
[4] C. Deffayet, Phys. Lett. B502, 199 (2001), hep-th/0010186. 



8 



P. J. E. Peebles and B. Ratra, Astrophys. J. 325, L17 (1988). 
R. R. Caldwell, Pliys. Lett. B545, 23 (2002), astro-ph/9908168. 

M. P. Dabrowski, T. Stacliowiak, and M. Szydlowski, Phys. Rev. D68, 103519 (2003), hep-th/0307128. 
A. Y. Kamenshchik, U. Moschella, and V. Pasquier, Phys. Lett. B511, 265 (2001), gr-qc/0103004. 

E. Witten (2000), hep-ph/0002297. 

D. J. C. MacKay, Information Theory, Inference and Learning Algorithms (Cambridge University Press, Cambridge, 2003). 
H. Akaike, IEEE Trans. Auto. Control 19, 716 (1974). 

H. Jeffreys, Theory of Probability (Oxford University Press, Oxford, 1961), 3rd ed. 
R. Trotta, Mon. Not. Roy. Astron. Soc. 378, 72 (2007), astro-ph/0504022. 

A. R. Liddle, P. Mukherjee, D. Parkinson, and Y. Wang, Phys. Rev. D74, 123506 (2006), astro-ph/0610126. 
G. Schwarz, Annals of Statistics 6, 461 (1987). 

W. Godlowski and M. Szydlowski, Phys. Lett. B623, 10 (2005), astro-ph/0507322. 
M. Szydlowski and W. Godlowski, Phys. Lett. B633, 427 (2006), astro-ph/0509415. 
R. Trotta, Mon. Not. Roy. Astron. Soc. 378, 819 (2007), astro-ph/0703063. 
M. V. John, Astrophys. J. 630, 667 (2005), astro-ph/0506284. 
M. Biesiada, JCAP 0702, 003 (2007), astro-ph/0701721. 

M. P. Hobson and C. McLachlan, Mon. Not. Roy. Astron. Soc. 338, 765 (2003), astro-ph/0204457. 
A. R. Liddle, Mon. Not. Roy. Astron. Soc. 351, L49 (2004), astro-ph/0401198. 

D. Parkinson, S. Tsujikawa, B. A. Bassett, and L. Amendola, Phys. Rev. D71, 063524 (2005), astro-ph/0409071. 
P. Serra, A. Heavens, and A. Melchiorri, Mon. Not. Roy. Astron. Soc. 379, 169 (2007), astro-ph/0701338. 
M. Szydlowski, A. Kurek, and A. Krawiec, Phys. Lett. B642, 171 (2006), astro-ph/0604327. 
M. Kunz, R. Trotta, and D. Parkinson, Phys. Rev. D74, 023503 (2006), astro-ph/0602378. 
A. R. Liddle, Mon. Not. Roy. Astron. Soc. Lett. 377, L74 (2007), astro-ph/0701113. 

M. Beltran, J. Garciar-Belhdo, J. Lesgourgues, A. R. Liddle, and A. Slosar, Phys. Rev. D71, 063532 (2005), astro- 
ph/0501477. 

P. Mukherjee, D. Parkinson, and A. R. Liddle, Astrophys. J. 638, L51 (2006), astro-ph/0508461. 

P. Mukherjee, D. Parkinson, P. S. Corasaniti, A. R. Liddle, and M. Kunz, Mon. Not. Roy. Astron. Soc. 369, 1725 (2006), 

astro-ph/0512484. 

A. Niarchou, A. H. Jaffe, and L. Pogosian, Phys. Rev. D69, 063515 (2004), astro-ph/0308461. 

T. D. Saini, J. Weller, and S. L. Bridle, Mon. Not. Roy. Astron. Soc. 348, 603 (2004), astro-ph/0305526. 

M. V. John and J. V. Narlikar, Phys. Rev. D65, 043506 (2002), astro-ph/0111122. 

D. Parkinson, P. Mukherjee, and A. R. Liddle, Phys. Rev. D73, 123523 (2006), astro-ph/0605003. 
A. Kurek and M. Szydlowski (2007), astro-ph/0702484. 
G. Efstathiou (2007), arXiv:0712.1513 [astro-ph]. 

S. Chongchitnan and G. Efstathiou, Phys. Rev. D76, 043508 (2007), arXiv:0705.1955 [astro-ph]. 
M. Szydlowski, A. Krawiec, and W. Czaja, Phys. Rev. E72, 036221 (2005), astro-ph/0401293. 
D. L. Dowe, S. Gardner, and G. Oppy, Br. J. Philos. Sci. 58, 709 (2007). 
W. Zimdahl, Int. J. Mod. Phys. D14, 2319 (2005), gr-qc/0505056. 

G. Olivarcs, F. Atrio-Barandcla, and D. Pavon, Phys. Rev. D74, 043521 (2006), astro-ph/0607604. 

G. Olivares, F. Atrio-Barandcla, and D. Pavon (2007), arXiv:0706.3860 [astro-ph]. 

M. Le Delliou, O. Bertolami, and F. Gil Pedro, AIP Conf. Proc. 957, 421 (2007), arXiv:0709.2505 [astro-ph]. 
M. Szydlowski, T. Stachowiak, and R. Wojtak, Phys. Rev. D73, 063516 (2006), astro-ph/0511650. 
K. Preese and M. Lewis, Phys. Lett. B540, 1 (2002), astro-ph/0201229. 

W. Godlowski, M. Szydlowski, and A. Krawiec, Astrophys. J. 605, 599 (2004), astro-ph/0309569. 

F. E. M. Costa, J. S. Alcauiz, and J. M. F. Maia (2007), arXiv:0708.3800 [astro-ph]. 
A. G. Riess et al, Astrophys. J. 659, 98 (2007), astro-ph/0611572. 

W. M. Wood-Vasey et al. (ESSENCE), Astrophys. J. 666, 694 (2007), astro-ph/0701041. 

T. M. Davis et al., Astrophys. J. 666, 716 (2007), astro-ph/0701510. 

Y. Wang and P. Mukherjee, Astrophys. J. 650, 1 (2006), astro-ph/0604051. 

D. J. Eisenstein et al., Astrophys. J. 633, 560 (2005), astro-ph/0501171. 

J. Simon, L. Verde, and R. Jimenez, Phys. Rev. D71, 123001 (2005), astro-ph/0412269. 

L. Samushia and B. Ratra, Astrophys. J. 650, L5 (2006), astro-ph/0607301. 

H. Wei and S. N. Zhang, Phys. Lett. B644, 7 (2007), astro-ph/0609597. 
J. Skilling, http://www.inference.phy.cam.ac.uk/bayesys/. 

A. Lewis, http://cosmologist.info/cosmomc/. 

A. Lewis and S. Bridle, Phys. Rev. D66, 103511 (2002), astro-ph/0205436. 
D. Parkinson, P. Mukherjee, and A. Liddle (http://www.cosmonest.org/). 



